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Abstract 

We consider a consistent extension of the SIMP models with dark mesons by including a dark 
U(1)_D gauge symmetry. Dark matter density is determined by a thermal freeze-out of the 3 —>• 2 
self-annihilation process, thanks to the Wess-Zumino-Witten term. In the presence of a gauge 
kinetic mixing between the dark photon and the SM hypercharge gauge boson, dark mesons can 
undergo a sufficient scattering off the Standard Model particles and keep in kinetic equilibrium 
until freeze-out in this SIMP scenario. Taking the S\J{Nf)xS\J{Nf)/Sl]{Nf) flavor symmetry 
under the SU(Vc) confining group, we show how much complementary the SIMP constraints on 
the parameters of the dark photon are for current experimental searches for dark photon. 
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I. INTRODUCTION 


Various evidences of dark matter (DM) imply that fundamental particles and interactions 


need to be extended beyond the Standard Model (SM) |l|. One of the appealing suggestions 
is the thermal DM scenario, where the DM relic density is determined through the freeze-out 
of the DM number changing process. Weakly Interacting Massive Particle (WIMP) provides 
the most popular thermal DM scenario, in which the annihilation of a DM pair into a pair 
of SM particles is responsible for the freeze-out. Since WIMP mass is of order weak scale 
for the effective coupling of aes ~ 0(10“^), ‘WIMP miracle’ has been the mainstream for 
thermal DM studies, corroborating the expectation of hnding new physics at the weak scale 
in the solutions for gauge hierarchy problem. 

Another interesting proposal for thermal (pseudo-) scalar DM has been recently made 


under the name of Strongly Interacting Massive Particle (SIMP) M. explaining DM relic 
density through the freeze-out of 3 —)■ 2 self-annihilation. The DM self-interaction is mo¬ 
tivated W potential small-scale problems J^, although it is strongly constrained by bullet 


cluster 3 and simulations on halo shape |7(. As a result, the SIMP scenario predicts dark 
matter with dimensionless self-interacting coupling of order one and mass in the 0.1 — 1 GeV 
range, which has not been explored seriously so far. 

The SIMP scenario requires the interaction between dark and SM sectors in the form 
of scattering for dark matter to be in kinetic equilibrium with the SM particles, without 
altering the structure formation 8|. Since such an inter-sector interaction also leads to the 
DM annihilation into SM particles, the inter-sector interaction strength is bounded from 
above for the dominance of 3 —)■ 2 self-annihilation, if combined, resulting in n.DM(o''i^)ann < 
n^y^{av‘^)3^2 < nsM{o'v)scatt, at the freeze-out temperature. Taking other constraints from 
ground-based experiments into account in addition, we can make quite a concrete prediction 
on the parameters of a specihc SIMP model. 


s. 


In this article, we consider a SIMP model with dark mesons suggested in Ref. [311, where 
the 5-point interactions between dark mesons for 3 —?■ 2 annihilation come from the leading 
interactions of the Wess-Zumino-Witten (WZW) term 9|,[l^. From the model point of view, 
the WZW term is interesting because it encodes various aspects on dark sector, namely, 
the WZW term exists only for a specihc havor symmetry of light dark quarks, depending 
on its spontaneous breaking pattern [ll], and it contains color number as a topological 
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index jJ. As an inter-sector interaction, we consider the hidden valley scenario 12|], that 
some heavy dark sector particle has a renormalizable coupling to the mediator particle that 
communicates between DM and SM particles. Higgs-portal interaction would be a natural 
candidate, but it is not enough for a sufficiently large DM-SM particle scattering at freeze- 
out temperature due to small Yukawa couplings of the light SM fermions. Thus, we study 
the case with a gauge kinetic mixing [^, that is the renormalizable and gauge invariant 
interaction between SM hypercharge U(l)y and dark sector U(1 )d with a dimensionless 
coupling. When the dark U(l) gauge coupling is not too tiny, dark meson annihilation into 
dark photons could easily dominate the annihilation process of dark matter. But, we can 
forbid it by taking the dark photon to be heavier than dark mesons. In this case, the gauge 
kinetic mixing plays a role of ‘hidden valley’ in the SIMP scenario. 

In Sec. m we briefly review dark mesons in the SIMP scenario, the abundance of which 
is frozen out by the WZW term. In Sec. IIIIl we discuss properties of the dark U(l) gauge 
symmetry, U(l)£), that is compatible with both the WZW term and the SIMP scenario. 
By taking SU(Yc) as an example of conhning gauge group, in Sec. HVl we present a viable 
parameter range for the dark photon mass and the strength of the gauge kinetic mixing. 
Finally, conclusions are drawn. 


II. SIMP DARK MESONS AND THE WZW TERM 


Dark meson appears as a composite state of dark quarks in models with dark conhning 


gauge group Gc and it has several interesting properties 


14| . The lightest mesons are 


interpreted as pseudo-Goldstone bosons from a spontaneous breaking of (accidental) havor 
symmetry, which guarantees their stability. Whereas havor symmetry is broken by higher 
dimensional operators, due to the compositeness of dark mesons, their decays induced by 
higher dimensional operator is suppressed, as compared to those of a fundamental scalar, so 
the stability of dark mesons is easier to achieve. For instance, dark mesons can be unstable 
by the decay into lepton pairs tt^ ii, due to the dimension-6 four-Fermi interaction. 




( 1 ) 


where M is the scale at which an explicit breaking of havor symmetry occurs, T“ is the 
generator of havor symmetry, and F is the dark meson decay constant of order of the con- 
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fining scale A. Note that this operator is diniension-6, rather than dimension-5, which is the 
dimension of the corresponding operator for a fundamental pseudo-scalar DM. Then, the 
lifetime of dark mesons is given by T”^ ~ 87rM^/(F^m^m^), and, for several hundred MeV 
to GeV-scale dark mesons and a similar conhning scale, tTo — )■ /i/I provides the largest decay 
rate. It is longer than the age of universe, 6.6 x lO^^GeV”^, as long as M is larger than 
10® GeV. Dark mesons can also decay into a pair of photons due to a dimension-7 operator, 
~ but the lifetime estimated is T"^ ~ 7rF^M®/(m^A®), 

resulting in a less stringent constraint on the cutoff, M > 10^ GeV. Moreover, the interac¬ 
tions between dark mesons induce the DM self-scattering, which provides a solution to the 
small-scale problems such as ‘core-cusp’ or ‘too-big-to-faih problems. 

In the SIMP scenario proposed in Ref. j^, the DM relic density can be explained from 
the freeze-out of dark mesons in the presence of their 5-point self-interactions, provided by 
the WZW term [s]. 




tNr, 


WZW 


2407r2 

N. I 


I dE^^’^^'^TrlU-^diUU-^djUU-^dkUU-^diUU-^dmU] 


( 2 ) 


2407r2 _ 

The WZW term is the outcome of a specihc flavor symmetry Gj and its spontaneous breaking 
to the subgroup H for a given conhning gauge group, relying on a nontrivial hfth homotopy 
group, 7r^{Gf/H) = 1 ll|. The manifest non-chiral global symmetry H is unbroken if it 
is respected by dark quark masses 15|. As a consequence, degenerate dark quark masses 
mg, or degenerate dark meson mass m^ = 8{A^/F‘^)mq guarantee the existence of the WZW 
term. 

Here, we quote the results obtained in Ref. j^, which will be useful for the discussion 
hereafter. The 3 —)■ 2 annihilation cross section is calculated from the WZW term to be 






(3) 


27r5Fio m \m 

where Tp is the freeze-out temperature, is the number of dark mesons, or dim{Gf/H), 
and F is a factor determined by group theory, ~ for large Nf. As a result of freeze- 
out, the number density of dark matter is given by udm = cTeqs/m.„, where Tgq = 0.8 eV 
is the matter-radiation equality temperature, s = (27r^/45)(/*5(T)T^ is the entropy den¬ 
sity of relativistic particles in thermal equilibrium, and c ~ 0.54 is the numerical con- 
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N„ 

*2 

N'ja? 

SU(Afc) 

SV{Nf)-X-SV{Nf) 

SU(iVj) 

(Nf > 3) 

Nj-1 

1 

1 

8{Nf - l){Nf + l){3Nj - 2Nj + 6) 

SO(Afc) 

S\]{Nf)/SO{Nf) 
{Nf > 3) 

^{Nf+2){Nf-l) 

^Nf{N]-l){N]-4) 

{Nf - l){Nf + 2){3Nj + 7Nj - 2Nj - 12Nf + 24) 

Sp(Afc) 

SV{2Nf}/Sp{2Nf) 
{Nf > 2) 

{2Nf + l){Nf - 1) 

lNf{N-j - 1){AN) - 1) 

4:{Nf - l){2Nf + l){6Nj - 7Nj - Nj + 3Nf + 3) 


TAB 

Ref. 


E I: Summary of group theory factors in the cases with nonzero WZW terms, quoted from 


a 


and Ref. 


111 . 


stant. For (crw^) 3^2 = freeze-out condition for the 3 —?■ 2 annihilation, 

n^y^{TF){(Tv‘^)^^ 2 {TF) = H{Tf)-, with the freeze-out temperature at Tp ~ mjr/20 { 2 ], deter¬ 
mines dark matter mass in terms of the effective DM self-coupling. 


— 0.03 Q!eff(TgqMp)^/^. 


(4) 


Thus, for Q!efr = 1 — 10, we get moM = 35 — 350 MeV. As will be discussed in a later 
section, in order to keep dark matter in kinetic equilibrium with heat bath, it is necessary 
to introduce the inter-sector interaction between dark and SM sectors. 

On the other hand, the leading 2^2 self-scattering comes from the kinetic term 
(F^/16)Tr(9^17(9^17“^), whose cross section is given by 

327rF4A'2’ ^ > 

where a? is another group theory factor ~ for large Nf. The self-interaction cross section 

is constrained to be CTseif/niTr < Icm^/g. This condition, together with the perturbativity 
bound of chiral perturbation theory, m-,^/F < 27r, imposes the dark meson mass to be in the 
0.1 — 1 GeV range, depending on the conhning gauge group. The group theory factors for 
possible gauge and flavor symmetries with nonzero WZW terms are summarized in Table [B 


III. DARK U(l) FOR SIMP DARK MESONS 

The dark U(l)p) charges of dark quarks are closely related to the form of quark mass 
terms. For SU(Ac) gauge groups, quarks and anti-quarks belong to fundamental and anti¬ 
fundamental representations, respectively, i. e. complex representations, so Dirac mass terms 
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such as {jnq)ijq^qj are allowed. In this case, dark quarks can be vector-like under U(1 )d so the 
model is automatically free from gauge anomalies. But, if U(1 )d is unbroken, dark mesons 
can be unstable in general, because dark mesons can decay fast into a pair of massless 
dark photons 7 ^ 5 , in the presence of AVV chiral anomalies. Even if the dark meson decays 
from AVV anomalies can be forbidden by appropriate 11(1)/) charge assignments, such as 
universal charges up to sign |^, we cannot prohibit a dark meson self-annihilation in the 
form of TTTT njo through AAAV anomalies^ 17|. Then, for the 3 —)■ 2 annihilation to be 
a dominant process for freeze-out, the dark gauge coupling for a unbroken U{1)d must be 
extremely small so the gauge kinetic mixing does not give an enough scattering cross section 
of dark mesons off the SM particles at freeze-out. 

For our later discussion on SU{Nc) confining groups, we take the U(l)/) compatible with 
SIMP dark mesons to be spontaneously broken so that dark photon gets massive. For dark 
photon mass my > rn.,r, the tttt Try/) processes from AAAV anomalies are kinematically 
forbidden. As will be discussed in the next section, in the presence of a gauge kinetic mixing 
between dark photon and the SM U(l)y, the off-shell processes, tttt — )■ vryj, — )■ Tre^e"*", opens 
up but it turns out to be suppressed as compared to the 3 —?■ 2 processes. 

For SO{Nc) and Sp{Nc) gauge groups, on the other hand, quarks in the fundamental 
representation, belong to real and pseudo-real representations, respectively, so there is no 
distinction between quarks and anti-quarks. As a result, only the Majorana mass terms are 
allowed. Denoting Weyl spinor indices as a, (3, - ■ ■, gauge multiplet indices as r, s, • • •, and 
flavor indices as dark quark mass terms appear as 


m, 




Qr,iQas,j + h.c.. 


( 6 ) 


in which = nigS^^S'^^ for SO{Nc) and for Sp(A^"c) gauge group, 

where J = ia 2 (D I is an antisymmetric second rank tensor. In this case, dark quarks are 
only chiral under U(l)/). Then, only after the U(l)/) is broken spontaneously, dark quarks 
obtain masses so does dark photon. 

When dark quarks are chiral under U(1 )d, a special care is needed. Suppose that a 
chiral U(l)/) is spontaneously broken and has a gauge kinetic mixing with the SM U(l)y, 
as for the case with a vector-like f/(l)/). First of all, there should be no gauge anomalies. 


^ Effects of both AVV and AAAV anomalies are encoded in the gauged WZW term 


1C 


m. 
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such as Gc — Gc—U(l)i:i or U(l)i:i-U(l)i:)-U(1)D anomalies. Secondly, even for dark photon 
mass with my > m.,^/2, the AVV anomalies could induce the decay of dark mesons into 
SM particles through the kinetic mixing, such as vr ^ 1*d1*d (e+e“)(e+e“). Therefore, 

the AVV anomaly terms should be forbidden. Finally, Partially Conserved Axial Current 
(PCAC) can couple to dark photon linearly such as Fd^nV^ in general. This results in the 
decay of dark mesons into a pair of SM fermions, so it is dangerous as well. These challenges 
with chiral 11 ( 1 ) 1 ) can be overcome by considering appropriate 11(1)/) charge assignments, 
possibly calling for extra heavy dark quarks. In our work, however, we won’t discuss this 
interesting case. 

In order to make the discussion simple, we consider vector-like dark quarks under 11(1))), 
that is spontaneously broken and has a kinetic mixing with SM U(l)y, and restrict ourselves 
to SU(Vc) conhning gauge symmetry. The physical results are not so different for the SO(Vc) 
and Sp(A^"c) gauge groups that need the breakdown of a chiral U(l))) for nonzero quark 
masses, as far as heavy dark quarks do not have order one Yukawa couplings. 

Dark mesons are pseudo-Goldstone bosons resulting from a spontaneous breaking of 
Gf =Sl] {Nf) xSlJ {Nf) down to H=S\J{Nf). As a minimal choice for a nonzero WZW 
term, we take Nf = 3. Furthermore, for the absence of the AVV anomalies, the U(l)/) 
charge operator Qd must satisfy Tr^QjyXa) = 0 with Xa being Gell-Mann matrices. Thus, 
we choose the charge matrix to be 


Qd 


^10 0 
0-10 
^0 0-1 


\ 

/ 


(7) 


Then, dark mesons are written as an SU(3)-valued matrix, U{x) = exp(2i where 




A 

F 


V2 


TT~ 

k- 


TT^ 

_L^o , j_~o 


iC+ 

k^ 


K 


w 


( 8 ) 


Due to the absence of AVV anomalies, neutral dark mesons are protected from the decays, 
27 £). One-loop corrections just rescale the vertices of AVV anomalies [^. In 
our case, the stability of neutral mesons is guaranteed by charge assignments in Eq. ([^ 
thanks to non-renormalization of AVV anomalies at all orders. Then, the kinetic term 
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/lQ)Ti{DfJJD^U~^), where the covariant derivative is D^U = d^U + igD[QDi U]V^, with 
qd being the dark gauge coupling (or od = fi'i)/47r being the dark structure constant), 
provides the leading interactions between dark mesons and dark photon V)(, 

= - i2gD{d^K^k~ - k^d^k~ + - 7r+(9^^")l/^ 

(9) 

+ Agl{k+k- + jf+%-)Vf,V>^. 

A remark on the effect of dark photon couplings on the mass splitting is in order. The 
U(1 )d charge assignment that we take makes some of dark mesons charged under U(l)£), 
resulting in a dark meson mass splitting coming from a£)A^Tr(Qf/(5f/“^), that violates flavor 
symmetry explicitly. If the dark meson mass splitting is large enough, the only lightest dark 
meson remains eventually as a result of SU(A"c) or U(1 )d interactions. However, U(l)£) 
mass contribution is small for a perturbatively small an- Since the SIMP scenario works for 
> 4 j^, for «£) = l/47r, the dark photon contribution to the mass splitting is as small 
as Ami < cdAVF^ oidF‘^ 0(10 ‘^)ml, i.e. less than 10%. Therefore, the dark meson 
mass degeneracy is a good approximation. Henceforth, we assume the U(1 )d charges given 
in Eq. (j7]). 


IV. SIMP DARK MESONS WITH SU(Ac) CONFINING GROUP 


We consider a gauge kinetic mixing between the U(l)i:i gauge boson (1%) and the U(l)y 
gauge boson (H^), given by 

1 1 sin V 1 

/:u(i)o = (10) 


After diagonalizing gauge kinetic and mass terms by 




Cw 

~SwC( + t^Si^ 

—Svi/S(; — 




K 

= 

Sw 

CwC( 

CwS( 



( 11 ) 

. . 


0 

-A 

Cx 

£C 





where 


tan 2C 


m%sw sin 2% 
ml - m|(c| - 


( 12 ) 


three mass eigenstates (A^, A'^) are interpreted as photon, Z-boson, and dark photon, 

respectively, with the masses of the latter two being 










(13) 


= 


2 L 


m 


"i|(l + s^4) + —^ ± _|_ 


^ 2 


4 

—m\mlr 


We get ~ m| and ~ my in the x 0 limit. To these gange bosons, electromagnetic 
(EM) cnrrent Jem, nentral Z—cnrrent J^, and dark sector cnrrent J)) conple, as 


■^int — 


+4l\ 


( —CvECftx)'^EM + + SW^x^d-^Z + 


(cveSC^x)'^EM + (^C Swtx^()Jz 

C\ 


(14) 


The leading interaction between dark and SM sectors is the dark photon conpling to EM 
current with shifted charges by {cwC(^tx )The kinetic mixing parameter, = cwc^tx, 
and the dark photon mass, my, are constrained by various experiments. 

For the range 0.02 GeV < my < 10.2 GeV, the recent BaBar experiment shows that 
G ~ 6 X 10“^ from the observation of e’''e“ 77 d* —^ 20|. The LHG experiments 

provide bounds, < 5 x 10“^ for 20 GeV < my < 30 GeV, from the analysis of a new Higgs 


decay mode, h —)■ Z'-jd, based on GMS8 21|, as well as e..^ < 10 ^ for 30 GeV < my < 70 GeV 
and my > 10^ GeV from the Drell-Yan (DY) ■jd production giving di-lepton signal 2^ [23 1 


based on Refs. [ 24 . These constraints are, however, imposed under the assumption that 
dark photon decays into SM particles only [2^. In our case, dark photon can decays mainly 


into a pair of dark mesons, if kinematically allowed, i.e. my > 2m.n-, so the decay branching 
fraction into a pair of visible SM particles, SMj, is modihed to be 

r(7,,^SMi) 


BR(7b ^ SM,) = 


r(7,, ^ SM,) 


- ■ 


(15) 


E, r( 7 B ^ SM,) + r(7z5 ^ tttt) r( 7 D ^ SM,) ’ 

where ^ SM,)/ ^ SM,) in the last equality corresponds to the old branching 

ratio without invisible decays. Since the experimental limits on the visible modes depend 
on e^Br( 7 E, e’''e“), the bound on gets weaker by a factor [aD/. Moreover, 
electroweak precision test (EWPT) provides stringent bounds, <2 x 10“^ for 10 GeV < 
my < 80 GeV, and e^y Y. 2.5 x 10“^ for my ~ m^, which is not affected by the invisible 
decays discussed above [23 1. 

Dark photon also couples to Z—current. For my <C m^, however, the Z—current coupling 
does not give a signihcant contribution, since the mixing angle approximated by C ~ —swX 
makes the coefficient for dark photon coupling to Z—current, ez = S(^ + swtx^C^ vanish at 
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the leading order. On the other hand, for the dark photon mass being around the Z—boson 
mass, the mixing angle gets larger as C ~ {rn\tw^'y )/cannot ignore it any 
longer and the U(l)o gauge boson is interpreted as a ‘dark Z boson’ 2^. This makes the 
lower bound on less stringent. 

We are now in a position to consider the conditions on my and that are consistent with 
the SIMP mechanism. At freeze-out temperature, Tp ~ ~ (5 — 50)MeV, photon, 

electron/positron, and neutrinos are relativistic particles in thermal equilibrium, and muon 
and pion begin to be non-relativistic. At that moment, the 3^2 self-annihilation from 
the WZW term is dominant over the other possible annihilation processes whereas dark 
meson-SM particle scattering processes do not decouple yet, provided that 


^DM(*^^)ann ^ )3—>-2 ^ ^SM(*^^)scatt) 


(16) 


at the freeze-out temperature, where the number density of a bosonic or fermionic SM 
particle is given by 


X 


27r2 ^ Jq ^\/x^+{msM/TF)'^ ^ ^ 

In order to prevent a pair of dark mesons from annihilating into a pair of dark photons 
due to the gauge interactions with we require my > m-,^. In this region, the 

TTTT —)■ SM SM annihilation rate, estimated as noM x [0{lQ‘^)aaD^‘^'rn\/{NT^my)], is smaller 
than the 3—7-2 annihilation rate, if 


< n ( '^v \ /0AG^\5/2 


(17) 


For typical parameters such as Tp ~ m7r/20, F ~ m^/5 and ttp ~ l/47r, the above condition 
is fulhlled for any value of satisfying the upper bounds given by ground-based experiments. 
The only exception appears around my = 2m^, where tttt —)■ SMSM annihilation rate is 
improved due to a resonance from l/(4m^ — my^. 

As for the dark meson scattering off the SM particles, the dominant process is tt -|- —?■ 
TT+e^ through the f—channel process, whereas ti+'J ti+'J is suppressed by a double kinetic 
mixing. Ignoring the lepton masses, we obtain the scattering cross section ioi tt + ^ ^ tt + ^ 
averaged over the number of dark mesons as 


((J'y)scatt,£ — 


4 

iW 


9 24(8s^ — 4sL -|- liei 

192e2 + ^ ^ ^ ^ ^ 


mz 


Tiaap 


m 


V 


T. 


m^ 


(18) 
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where a factor 4 represents the degrees of freedom of U(l)i:i charged mesons, and tt^. On 
the other hand, the scattering cross sections of dark mesons off the neutrinos, vr +—)■ tt + z/ 
and the SM pion, vr + tt^j^ ^ tt + are given by 


{av) 


scatt.z/ 


4 247raa£)e^?Ti^ (Tp 


{av) 


SCatt.TT 




, 192el 1, _ 

+ ;:rT#7(l - 24) 


ml fTp 


TTO'O' p) —j 


rriy 


From the condition. 


(19) 


i=e,fi 


scatt,7r5 


( 20 ) 


the larger my, the stronger the lower bound on gets, according to 




OtD^"! I / 
Vmy / 


( 21 ) 


We present constraints on my and for several cases with different conhning groups 
in Fig. [1] and Fig. O where ap = l/dvr and an = 0.01 were taken, respectively, and 
that dark meson masses are assumed to be degenerate. For W = 4, the minimal W that 
the SIMP mechanism works, only m^ ~ 0.45 GeV is allowed because the perturbativity 
condition, X = m-^/F < 2%, and the self-interaction bound, aseii/m-„ < Icm^/g, almost 
coincide |3|. For W = 6 and W = 10, a wider range of dark meson masses are allowed 
such as 0.37 GeV < m^r < 0.56 GeV and 0.26 GeV < m^ < 0.8 GeV, respectively, and x is 
hxed by the DM relic density. In both cases, the upper bounds satisfy x = 2 ti and the lower 
bounds satisfy x = 5.48 and 4.6, respectively. 

We note that dark photon is taken to be heavier than dark mesons in order for the 3—7-2 
annihilation to dominate over the tttt 'Jd'Id annihilation. As a result, there appears a 
lower bound, e.y > 10“^ at my = m^r, due to Eq. fl2ll) . On the other hand, the AAAV 
anomalies induce the annihilation process such as 7i{ki)7i{k2) '^{k 3 )e~{pi)e'^{p 2 ) through 
off-shell dark photon, where momenta of particles in the process are explicitly written. The 
interaction vertex contains e^ppakik 2 k^v{p 2 )'y^u{pi) term, which vanishes in non-relativistic 
limit, /ci ~ /c 2 — (m7r,0). Thus, the annihilation cross section for this process is estimated 
as aa£)el^m^Tp/{NT,myF^). Therefore, as compared to the annihilation cross section of 
TTTT —)■ which is smaller than the one for the 3 —)■ 2 processes, the annihilation cross 
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Logjo[mv (GeV)] 


Logjo[mv (GeV)] 


FIG. 1: Bounds on my vs for dark mesons being compatible with the SIMP scenario for I/dvr. 
Three figures correspond to Gc =SU(4), SU(6), and SU(IO), respectively. Imposed constraints, 
distinguishable by colors, are written explicitly, while the allowed parameter space is uncolored. 
For my > 2m,r) BaBar and LHC bounds are rescaled taking —)• 27r invisible decay into account. 


section of vrvr —)■ 7re“e+ is suppressed by [(27r)^F®] ~ m®/(27rF)®, within the valid 

regime of chiral perturbation theory. 

The lower bound on for a given my follows from the estimation given by Eq. (j2Tl) . and 
the upper bound comes from ground-based experiments. We also find from Eq. fl2T]) that 
the SIMP condition requires a large kinetic mixing for a large dark photon mass, eventually 
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Logjo[my (GeV)] 



Logjo[mv (GeV)] 


Logjo[mv (GeV)] 


FIG. 2: Bounds on my vs for dark mesons being compatible with the SIMP scenario for 

Oi£) = 0 . 01 . 


constrained by ground-based experiments. In summary, for N^. < 10, dark photon masses 
are allowed up to my ~ 10^ GeV with varying limits on e.^. We note that there are more 
allowed values of e.^ around my ~ mz due to the non-negligible contribution from ez- The 
bounds get stronger near my = 2m^, where the dark meson annihilation into a pair of SM 
particles becomes enhanced as discussed previously. 
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V. CONCLUSION 


We have considered an extension of the models with SIMP dark mesons by including a 
dark local U(l)£) symmetry under which dark quarks are vector-like. Dark mesons are still 
good candidates for SIMP DM, as the chiral anomalies associated with U(l)£) are absent. In 
the presence of a gauge kinetic mixing between U(l)£) and the SM U(l)y, the dark sector is 
communicated with the SM particles through the Z' portal so that it can be kept in kinetic 
equilibrium with the SM sector until the freeze-out in the SIMP scenario. 

The SIMP conditions restrict the parameter space for dark photon mass my and kinetic 
mixing that is otherwise unconstrained by ground-based experiments. Focusing on dark 
mesons living on the SU{Nf)xSU{Nf)/SU{Nf) flavor symmetry and taking the SU(iVc) 
conhning group for them, we showed that the combination of the SIMP conditions with 
various ground-based experiments searching for dark photon can restrict the parameter 
space to < my < 10^ GeV and 10“^ < e.y < (10“^ — 10“^), for dark gauge coupling of 
order one and W < 10. 

Acknowledgment s 

We would like to thank Su Min Choi, Ji-Haeng Huh, Kwang Sik Jeong, Hye-Sung Lee 
and Jong-Chul Park for discussions. The work of HML is supported in part by Basic Science 
Research Program through the National Research Foundation of Korea (NRF) funded by 
the Ministry of Education, Science and Technology (2013R1A1A2007919). The work of MS 
is supported by IBS-R018-D1. 


[1] G. Bertone, D. Hooper and J. Silk, Phys. Kept. 405, 279 (2005) [hep-ph/0404175]. 

[2] Y. Hochberg, E. Kuflik, T. Volansky and J. G. Wacker, Phys. Rev. Lett. 113, 171301 (2014) 
[arXiv: 1402.5143 [hep-ph]]. 

[3] Y. Hochberg, E. Kuflik, H. Murayama, T. Volansky and J. G. Wacker, arXiv:1411.3727 [hep- 
ph]. 

[4] N. Bernal, C. Garcia-Gely and R. Rosenfeld, arXiv:1501.01973 [hep-ph]. 


14 



[5] D. N. Spergel and P. J. Steinhardt, Phys. Rev. Lett. 84, 3760 (2000) [astro-ph/9909386]; 
W. J. G. de Blok, Adv. Astron. 2010, 789293 (2010) [arXiv:0910.3538 [astro-ph.CO]]; 
M. Boylan-Kolchin, J. S. Bullock and M. Kaplinghat, Mon. Not. Roy. Astron. Soc. 415, 
L40 (2011) [arXiv:1103.0007 [astro-ph.CO]]; D. H. Weinberg, J. S. Bullock, F. Governato, 
R. K. de Naray and A. H. G. Peter, arXiv: 1306.0913 [astro-ph.CO]. 

[6] M. Markevitch, A. H. Gonzalez, D. Clowe, A. Vikhlinin, L. David, W. Forman, C. Jones 
and S. Murray et ai, Astrophys. J. 606, 819 (2004) [astro-ph/0309303]; D. Clowe, A. Gon¬ 
zalez and M. Markevitch, Astrophys. J. 604, 596 (2004) [astro-ph/0312273]; S. W. Randall, 
M. Markevitch, D. Clowe, A. H. Gonzalez and M. Bradac, Astrophys. J. 679 (2008) 1173 
[arXiv:0704.0261 [astro-ph]]. 

[7] M. Rocha, A. H. G. Peter, J. S. Bullock, M. Kaplinghat, S. Garrison-Kimmel, J. Onorbe and 
L. A. Moustakas, Mon. Not. Roy. Astron. Soc. 430, 81 (2013) [arXiv:1208.3025 [astro-ph.CO]]; 
A. H. G. Peter, M. Rocha, J. S. Bullock and M. Kaplinghat, Mon. Not. Roy. Astron. Soc. 
430, 105 (2013) [arXiv:1208.3026 [astro-ph.CO]]. 

[8] E. D. Carlson, M. E. Machacek and L. J. Hall, Astrophys. J. 398, 43 (1992); A. A. de Laix, 
R. J. Scherrer and R. K. Schaefer, Astrophys. J. 452 (1995) 495 [astro-ph/9502087]. 

[9] J. Wess and B. Zumino, Phys. Lett. B 37, 95 (1971). 

[10] E. Witten, Nucl. Phys. B 223, 422 (1983). 

[11] E. Witten, Nucl. Phys. B 223, 433 (1983). 

[12] M. J. Strassler and K. M. Zurek, Phys. Lett. B 651, 374 (2007) [hep-ph/0604261]. 

[13] L. B. Okun, Sov. Phys. JETP 56, 502 (1982) [Zh. Eksp. Teor. Eiz. 83, 892 (1982)]; B. Holdom, 
Phys. Lett. B 166, 196 (1986). 

[14] J. M. Cline, Z. Liu, G. Moore and W. Xue, Phys. Rev. D 90, no. 1, 015023 (2014) 
[arXiv:1312.3325 [hep-ph]]. 

[15] C. Vafa and E. Witten, Nucl. Phys. B 234, 173 (1984); C. Vafa and E. Witten, Commun. 
Math. Phys. 95, 257 (1984). 

[16] O. Bar and U. J. Wiese, Nucl. Phys. B 609, 225 (2001) [hep-ph/0105258]. 

[17] W. A. Bardeen, Phys. Rev. 184, 1848 (1969). 

[18] O. Kaymakcalan, S. Rajeev and J. Schechter, Phys. Rev. D 30, 594 (1984); N. K. Pak and 
P. Rossi, Nucl. Phys. B 250, 279 (1985). 

[19] J. F. Donoghue and D. Wyler, Nucl. Phys. B 316, 289 (1989); B. Borasoy and R. Nissler, 


15 


Eur. Phys. J. A 19, 367 (2004) [hep-ph/0309011]. 

[20] J. P. Lees et al. [BaBar Collaboration], Phys. Rev. Lett. 113, no. 20, 201801 (2014) 
[arXiv: 1406.2980 [hep-ex]]. 

[21] D. Curtin, R. Essig, S. Gori, P. Jaiswal, A. Katz, T. Liu, Z. Liu and D. McKeen et ai, Phys. 
Rev. D 90, no. 7, 075004 (2014) [arXiv: 1312.4992 [hep-ph]]. 

[22] J. M. Cline, G. Dupuis, Z. Liu and W. Xue, JHEP 1408, 131 (2014) [arXiv: 1405.7691 [hep- 
ph]]; 1. Hoenig, G. Samach and D. Tucker-Smith, Phys. Rev. D 90, no. 7, 075016 (2014) 
[arXiv: 1408.1075 [hep-ph]]. 

[23] D. Curtin, R. Essig, S. Gori and J. Shelton, JHEP 1502, 157 (2015) [arXiv:1412.0018 [hep- 
ph]]. 

[24] S. Ghatrchyan et al. [CMS Collaboration], JHEP 1312 (2013) 030 [arXiv:1310.7291 [hep-ex]]; 
[ATLAS Collaboration], ATLAS-CONP-2013-017, ATLAS-COM-CONP-2013-010. 

[25] B. Batell, M. Pospelov and A. Ritz, Phys. Rev. D 79, 115008 (2009) [arXiv:0903.0363 [hep- 
ph]]. 

[26] H. Davoudiasl, H. S. Lee and W. J. Marciano, Phys. Rev. D 85, 115019 (2012) [arXiv: 1203.2947 
[hep-ph]]. 


16 


